Phonon-phonon interactions and phonon damping in carbon nanotubes 
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We formulate and study the effective low-energy quantum theory of interacting long-wavelength 
acoustic phonons in carbon nanotubes within the framework of continuum elasticity theory. A gen- 
eral and analytical derivation of all three- and four-phonon processes is provided, and the relevant 
coupling constants are determined in terms of few elastic coefficients. Due to the low dimensional- 
ity and the parabolic dispersion, the finite-temperature density of noninteracting flexural phonons 
diverges, and a nonperturbative approach to their interactions is necessary. Within a mean-field 
description, we find that a dynamical gap opens. In practice, this gap is thermally smeared, but 
still has important consequences. Using our theory, we compute the decay rates of acoustic phonons 
due to phonon-phonon and electron-phonon interactions, implying upper bounds for their quality 
factor. 

PACS numbers: 63.22.Gh, 62.25.-g, 63.20.kg, 63.20.kd 
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INTRODUCTION 



Even after more than a decade of very intense re- 
search efforts,— the unique electronic and mechanical 
properties of carbon nanotubes (CNTs) continue to at- 
tract considerable interest. A major driving force for 
this interest comes from the prominent role played by 
phonons in CNTs. Phonons are crucial when interpret- 
ing experimental data for resonant Raman or photo- 
luminescence excitation spectra^i^ and for the under- 
standing of electrical^ and thermal^ transport in CNTs. 
Moreover, phonons are responsible for interesting na- 
noelectromechanical effects in suspended CNTsj^iii^iiSiii^ 
and they lead to quantum size effects in the specific 
heatjii The real-time nonlinear dynamics of a CNT 
phonon mode has also been monitored experimentally by 
femtosecond pump-probe techniques (coherent phonon 
spectroscopy) J^ii^ 

Recent experiments have shown that mechanical os- 
cillations of suspended carbon nanotubes can be ex- 
cited by a cantilever and detected by scanning force 
microscopy] ^ ^ Such experiments yield both the fre- 
quency oj and the quality factor Q = lo /V (with de- 
cay rate V) of the respective phonon mode. The so 
far observedi^ values, Q < 10"*, imply significant de- 
cay rates even at rather low temperatures, and require 
to identify the relevant decay channels for phonons in 
individual CNTs. Our paper is primarily devoted to un- 
derstanding the importance of phonon-phonon interac- 
tions in such decay processes. The quality factor can also 
be extracted from Raman spectroscopy^ and from nano- 
electromechanical measurements, using phonon-assisted 
Coulomb blockade spectroscopy^ or capacitive detec- 
tion of mechanical oscillations4 In principle, coherent 
phonon spectroscopy also allows to access damping rates 
of phonon modes, and hence their quality factors. Very 
recently, the possibility of cooling a vibrating carbon 
nanotube to its phononic ground state has also been 



discussedfi^ 

The recent experimental progress described above 
highlights the need for a reliable theory of phonon- 
phonon (ph-ph) interactions in CNTs. On the theoret- 
ical side, many authors have analyzed the noninteract- 
ing problem, i.e. the harmonic (or linear) theory, which 
allows to derive explicit theoretical results for the ther- 
mal conductanceiii^ and for the specific heatJ^ Moti- 
vated by the observation that molecular dynamics cal- 
culations seem to be in good agreement with thin-shell 
model predictions^^ several theoretical works2ii22ii^i2ii^ 
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have adapted thin-shell hollow cylinder model: 
the calculation of phonon spectra. However, the thin- 
shell approach leaves open the question of how to actu- 
ally choose the width of the carbon sheet. A popular 
and more microscopic approach is to instead start from 
force-constant modelS )^^i^^'^° taking into account up to 
fourth-nearest-neighbor couplings in the most advanced 
formulations j ^ ^ i These calculations predict four acous- 
tic phonon branches (with uj{k 0) = 0), namely a 
longitudinal stretch mode, a twist mode, and two degen- 
erate flexural modes (see Sec. Ill CI for details). The re- 
sulting phonon spectra are in very good agreement with 
a much simpler calculation based on continuum elastic- 
ity theoryf^^i^ building on the known elastic isotropy of 
the honeycomb latticei^ The elasti c ap proach will be 
employed in our study as well. Ref. |36| provides a gen- 
eral discussion of the accuracy of elastic continuum the- 
ories for phonons in CNTs. For very small CNT radius 
R, however, hybridization effects involving carbon a or- 
bitals lead to qualitative changes, elastic continuum the- 
ories (at least in the form below) may break down, and 
first-principle calculations become necessary<^2i2^ 

In contrast, the problem of ph-ph interactions in CNTs 
is much more difficult and has been treated in only a 
few works, although phonon anharmonicities are impor- 
tant for several physical observables^^i^i^ e.g. for ther- 
mal expansion (which is a controversial issue in CNT 
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theorji^), to explain the stability of low-dimensional ma- 
terials (which would be unstable in harmonic approxi- 
mation due to the Mermin- Wagner theorem), in order 
to establish a finite thermal conductivity, or to pro- 
vide a finite lifetime for phonons. The latter issue is 
particularly relevant in our context, but aside from a 
numerical high-temperature study^ which ignored the 
(lowest-lying) fiexural phonons, to the best of our knowl- 
edge ph-ph interactions in CNTs have only been studied 
by Mingo and BroidoJ^iii Their work considered three- 
phonon processes and their effects within a Boltzmann 
transport equation. The main conclusion of Refs. ITsl . 
H was that anharmonic effects are generally weak but 
important in establishing upper bounds for the thermal 
conductance. In addition, they computed the lengthscale 
up to which phonons show ballistic motion. Where ap- 
plicable, our results below are in accordance with theirs, 
but four-phonon processes (which govern the decay of 
flexural phonons) have not been studied so far. 

We shall consider two important mechanisms for the 
decay of long-wavelength acoustic phonons in single- wall 
CNTs, namely electron-phonon (e-ph) and ph-ph scatter- 
ing. We show that the dominant e-ph coupling terms (re- 
sulting from the deformation potential contribution) do 
not allow for phonon decay due to kinematic restrictions, 
and thus an intrinsic upper bound for the temperature- 
dependent quality factor of the various modes can be 
derived from ph-ph interactions alone. These upper 
bounds are given below. The problem of phonon de- 
cay has in fact a rather long history. Early work on the 
decay of an optical phonon into two acoustic phonons 
via anharmonicitie a^^i^^ proposed a scheme for nonlin- 
ear phonon generation. Phonon decay via ph-ph inter- 
action is also important for the understanding of neu- 
tron scattering data^ and for the collective excitations 
in liquid helium4i Such effects have even been consid- 
ered in a proposal for a phonon-based detector of dark 
matter.— General kinematic restrictions often prevent 
the decay of phonon modes. Lax ei aZ. have shown^ that 
a given acoustic phonon cannot decay into other modes 
with higher velocity at any order in the anharmonicity. 
For the lowest-lying acoustic phonon mode, one then ex- 
pects anomalously long lifetimes, while the higher acous- 
tic modes typically decay — in three-dimensional (3D) 
isotropic media with rate T cx |p|^ for phonon momen- 
tum p.i^i^ Such questions are particularly interesting 
in the CNT context, where a degenerate pair of flexural 
modes has the lowest energy, and the low dimensional- 
ity and the quadratic dispersion relation of the flexural 
mode may give rise to unconventional behavior. 

Before describing the organization of the paper, we 
pause for some guidance for focused readers. Experimen- 
tally minded readers can find our central predictions for 
the decay rate (and hence the quality factor) of the low- 
energy phonon modes in Eqs. ((54)) . (|58|) and ((67)) . The 
dependence of the resulting Q factors on the CNT radius 
R is shown in Fig. [21 Those interested in the main new 
theoretical results will find them in Eqs. ((^S)) and (pS)) . 



where the complete low-energy Hamiltonian for interact- 
ing acoustic phonons in single-wall CNTs is given, with 
the modified flexural dispersion relation (|l5)l . This mod- 
ification takes into account the instability of a harmonic 
theory implied by the Mermin- Wagner theorem, and in- 
cludes interaction effects in a nonperturbative manner. 
The calculation of the decay rates is then possible in 
a perturbative manner, and leads to the results quoted 
above. 

Let us conclude this Introduction with the organization 
of the paper. In this work, based on the elastic contin- 
uum description, we formulate a complete and analytical 
theory of interacting long-wavelength acoustic phonons 
in single- wall CNTs. In Section |TT] we show that the sim- 
plicity of the elastic approach allows us to go beyond 
the harmonic approximation (which is briefly reviewed 
in Sec. Ill C} . and thereby provides a complete theory of 
all possible three- and four-phonon scattering processes, 
described in detail in Sec. IIIII The theory is then ap- 
plied in Sec. [IV]to the calculation of phonon decay rates. 
We thereby infer intrinsic upper bounds for the quality 
factor of the relevant acoustic modes. We comment on 
effects of e-ph interactions on the quality factor in SecfVl 
and end the paper with a discussion and an outlook in 
Sec. IVH Calculational details have been relegated to two 
appendices. Finally, we note that while some of our re- 
sults are also relevant to 2D graphene monolayers ^^i*^!^ 
for the sake of clarity, we restrict ourselves to the CNT 
case throughout the paper. We sometimes set h ^ I in 
intermediate steps. 

II. NONLINEAR STRAIN TENSOR AND 
ELASTIC THEORY 

In this section we shall develop the low-energy the- 
ory of interacting long-wavelength acoustic phonons in 
CNTs. To be specific, we first discuss semiconducting 
single-wall CNTs, where e-ph scattering processes can 
safely be ignored. 

A. Strain tensor in cylindrical geometry 

We start from a continuum description, where 
long-wavelength phonons are encoded in the three- 
dimensional displacement field, u, with local-frame com- 
ponents Un=x,y,z (see below). The surface of an unde- 
formed cylinder, representing the CNT with radius R, is 
parametrized as 

'R.{v) = Re^{x) + yey, v^{x,y). (1) 

We use cylindrical coordinates with x/ R (where Q < x < 
2-kR) denoting the angular variable. The corresponding 
local- frame unit vector is Bx{x), while points along the 
cylinder axis and ez{x) is perpendicular to the cylinder 
surface, i.e. z corresponds to the radial coordinate. Note 
that Rdx^z = and Rdx^x — —e^, and R depends only 
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on the coordinates r = (x, y) parametrizing the cyHnder 
surface. Our convention for the coordinates follows the 
notation of Ref. [13, which is convenient because it con- 
nects the problem on the cylinder (CNT) to the one on 
the plane (graphene). 

The surface of the deformed cylinder is then 
parametrized in terms of the displacement field as 

x(r) = R(r) + u(r) = R(r) + ^ u„(r)e„(x). (2) 

n—x.y,z 

Equations ([1]) and ^ imply the relation 



dux 



R 



dx 



e^ix) 



(3) 



+ [duy + dy\ey + 
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dUz f^d.x 

R 



where contributions come both from the variation of the 
displacement field and from the change in the local frame. 
Given the displacement field, the symmetric strain tensor 
Uij{v), with i,j = X, y, can be obtained from the defining 
relation^^ 

Uxxdx^ + Uyydy^ + 2uxydxdy — '^{dy^ — dR^). (4) 

Employing Eq. ([3]), after some algebra, the strain tensor 
follows. It is composed of a linear and a nonlinear part 
in the displacement field, u = u^' 



2m: 



nlin 



DiUj 



{DiUn){D.jUn), 



(5) 
(6) 



where we use covariant derivatives, 



n dux Uz duz Ux 



while DyUn 



dx "-V 



g^M„ and DxUy 
One easily verifies that the strain tensor Uij respects 
fundamental symmetries. In particular, = for arbi- 
trary rigid translations or rotations of the whole cylinder. 
For translations, both the linear and the nonlinear part 
of the strain tensor vanish separately, but this is not the 
case for rotations. While u!'" = under infinitesimal ro- 
tations, the full nonlinear strain tensor must be kept in 
order to correctly account for Uij = under finite rota- 
tions. 



to u„. The theory is quantized via the standard commu- 
tation relations [with r = (x, y) and n, n' = x, y, z]. 



[p„(r),u„/(r')]_ = -ihSnn'S{r - r'). 



(9) 



Armed with the nonlinear strain tensor, progress is now 
possible by invoking symmetry considerations and the 
usual assumption of a space-time local elastic energy 
density U = U[u(r,t)] depending only on the strain ten- 
sor. The clastic energy density is then expanded in the 
strain tensor up to fourth order, U[u\ = U2 + U3 + U4. 
This expansion will fully account for all elementary ph- 
ph scattering processes involving at most four phonons. 
In order to give explicit expressions, we will exclude 
the case of ultrathin CNTs, which is difficult to model 
with an elastic continuum approach. For extremely 
small radius, a — n orbital hybridization effects due to 
the curvature of the cylinder can lead to dramatic ef- 
fects and, in particular, may change the honeycomb lat- 
tice structure.— In practice, this means that we require 
R > 4A. In that case, curvature effects generally scale 
as and as outlined in Appendix [Xl their inclu- 

sion is possible on phenomenological grounds within our 
nonlinear elasticity theory. Ignoring curvature effects for 
the moment, a straightforward connection to the corre- 
sponding planar problem of graphene can be established. 
Since graphenc's honeycomb lattice is isotropic with re- 
spect to elastic properties,— U[u] can only depend on 
invariants of the strain tensor under the symmetry group 
0(2). Independent invariants can then be formed using 
the trace of the strain tensor, Trw, and its determinant, 
det u = [(Tru)2 - Ttu^]/2. 

Starting with U2[u\, quadratic in the strain tensor, one 
arrives at the familiar expression^ 



-2 (T™^' 



/i Tru 



(10) 



where A and /i are Lame coefficients. Their value in 
graphene is estimated to be (see, e.g., Refs. ISSIIH^ ) 

K o m^ a a 
-2.90x10®^, — -1.51x10^^, (11) 

PM pM 

with the bulk modulus K = p + X. The 2D Poisson ratio 
then corresponds to 



K - n 
V = ~ 0.31, 



(12) 



B. Elastic energy density 

The Hamiltonian density is given by the sum of the 
kinetic and the elastic energy density, 



(8) 



where p^ = 3.80 x 10 ^ kg/m^ is the mass density of 
graphene, and p„ is the canonically conjugate momentum 



in agreement with that computed using an empirical 
force-constant modeljS Note that is already non- 

linear in the displacement field due to the nonlinearity 
^ of the strain tensor. We shall refer to such nonlinear- 
ities, resulting already fromZ//2, as geometric. Geometric 
nonlinearities do not involve new material parameters be- 
yond the Lame coefficients. 

Apart from geometric nonlinearities. there are also an- 
harmonic contributions (W3 and due to higher-order 
terms in the expansion of the elastic energy density in the 
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strain tensor. In cubic order, one can build three invari- 
ants from the strain tensor, namely Tru'^, Tru^ Tru, and 
{Tru)^. However, these invariants arc not independent, 
since 2 Trit^ = 3 Tru Tru^ — (Trit)^. Hence there are just 
two new anharmonic couplings in third order, denoted as 
^1 and ^2, leading to 

Usiu] = ^liTruf + ^2 Trw^ Tru. (13) 

Taking u = u^'", this produces three-phonon interaction 
processes from anharmonic terms in the elastic energy 
density, on top of the geometric nonlincaritics. Note 
that the nonlinear part u"''" then causes four-phonon 
processes (or higher orders) from Eq. (jl3|l . Finally, in 
quartic order there arc five invariants, 

Tiu^, Tru^ Trit, (TrM^)^, Tiu^ (Trw)^ (Tvu)'^. 

However, because of the identities 

2Trii'* = -{Ttu)'^ + {Tru^ f + 2{Ttu)^ Ttu^, 

Tru Tru^ = Tru"^ - i {(Ttu^)^ - {Truf Ttu^) , 

only three out of the five invariants are independent. 
With fourth-order anharmonic couplings (ki, K2, k^), wc 



can thus write 

lli[u] = Ki(Tru)^ + K2(Tru)2 Tru^ + KsiTiu^f. (14) 

As we show below, cf. Eqs. ([Ml), and ([671, the dom- 
inant decay processes for acoustic phonons are governed 
by the geometric nonlinearities alone, and no parameter 
estimates for the anharmonic couplings (ki,2,3 and ^1,2) 
are necessary for the calculation of these decay rates. 
This remarkable result could not have been anticipated 
without explicit computation of all contributions. 

The final step is to insert u = u''" + m"''" in U[u], 
and thereby to separate the harmonic theory (nonin- 
teracting phonons, Hq) from interactions (three-phonon, 
H''^\ and four-phonon, H^^\ processes), where H = 
Ho + H^"^^ + H'^^\ Wc do not take into account higher- 
order ph-ph scattering processes beyond the fourth order. 
Collecting terms, the harmonic theory corresponds to the 
Hamiltonian density 

+ \ (T™''")' + Tr [{u^n'] . (15) 

n 

All possible three-phonon processes are encoded in 



H^^^ = A Trw""" Tru''" + 2^l Tr(w""u""") + (Tru"") + ^2 Tru"" Tr [(u"")^] , 
while all four-phonon processes are contained in 

n^^^ = - (Tru""")^ ^ Tr [(u""")^] + 3^1 Tru""" (Tru"")^ + ^Tru"" [2Tr(u""u""") + Tru""" Tru""] 
+ Ki(Tru"")4 + K2(Tru"")2 Tr [(u"")^] + ^3 (Tr [(u"")^])' . 

I 



(16) 



(17) 



While this may seem like a rather complicated theory, 
we shall see below that the geometric nonlinearities (i.e. 
the terms corresponding to the Lame coefficients A and 
/i) already generate the most relevant structures. 



C. Harmonic theory 

Let us first diagonalize the noninteracting Hamiltonian 
Hq, see Eq. ([T5|) . and thereby determine the phonon spec- 
trum. Although the results of this subsection have essen- 
tially been obtained before,— 1^ we repeat the main steps 
in order to keep the paper self-contained. First, we per- 
form a Fourier transformation of the displacement field 
u„(r), introducing the momentum hk along the y-axis 
and the integer angular momentum quantum number ^, 



1 



V2^ 



^e*^^/-^+*'=^u„(fc,^). 



where u|,(fc,£) = u„(-fc,-€) and Y.k = ^nd 
an analogous transformation for p„. The commutation 
relations (O then read 

[pn{k,t),un'{-k',-i')y = -2mh5nn'5u,5{k-k'). 

Some algebra yields H^ in the form 

Ho^ f dxdyHo = J-Y^pl{k,e)p^ik,£) (18) 

+ I ui{kJ)Ann'{k,e)u„,{kJ), 



M 



where the elastic matrix A(fc,£) = (A„„')(A:,^) is given 

by 



A = 



/ I'iK+ti.) , 2 MK ii{K+p.) 

^ + (i^ + ^)fc2 ^^HK-^ 
. il(K+ti) ik{K-fj,) k+fj, 



k,£ 



(19) 
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This 3x3 matrix is obviously Hcrmitian and obeys the 
time-reversal symmetry relation^ A{—k, —£) = A*(fc, £), 
where the star denotes complex conjugation. Note that 
the chirality of the CNT does not affect the clastic matrix 
(and hence the dispersion relation) within the low-energy 
theory. However, the situation is different for high-energy 
optical phonons or when taking e-ph interactions into 
account.^ 

The normal- mode frequencies LUj(k,i) with corre- 
sponding polarization unit vectors ej{k,£) (the index J 
labels the normal modes) then follow from diagonalizing 
the elastic matrix, 



A(k,e)ej{k,e) = pMUjj{k,£)ej{k,£). 



(20) 



The above symmetries of the elastic matrix imply 
wj(-fc,-£) = wj(fc,^) and [ej{k,£)]* = ej{-k,-l). 
Moreover, polarization vectors for given k and £ are or- 
thonormal, e*j{k,£) ■ ej,{k,£) = 6jj'. Expanding the dis- 
placement field in terms of the polarization vectors and 
introducing boson creation, aj{k,£), and annihilation. 



j{k, £), operators, 

[aj{k,£),a''j,{k',f)]_ = 2Ti5{k - k')5.,j.5t 



(21) 



we arrive at the quantized noninteracting phonon Hamil- 
tonian, 



Ho = ^hwj{k,£) ia}j{k,£)aj{k 



Jki 



(22) 



The displacement field components arc then 
«„(r) = -^^e^^-/^+^'=^[e^(fc,£).e„]«,,(fc,^), (23) 



Jki 



with the normal-mode components^ expressed in terms of 
the boson operators, 



uj{k,£) = 



2pMUJj{k,£) 



(^aj{k,£)+al{-k,~£)^ 



(24) 

We next summarize the solutions of the eigenvalue prob- 
lem (PU)) . We are interested in the long- wavelength 
(|fc|i? ^ 1) phonon modes, in particular those with 
Luj{k -^0,£)^ 0. 

In the £ = sector, there arc three eigcnmodes, namely 
J ^ T (twist mode), J = L (longitudinal stretch mode), 
and J = B (breathing mode). For the twist mode, we 
find for arbitrary k the result 



erik) 




(25) 



where vt = 1-23 x lO^m/s. Note that e^; points along 
the circumferential direction. For the longitudinal stretch 



mode, we obtain 



^L{k) = WL|fc|-f 0(fc2), VL 



4Kn 



Pm{K + p) 



(26) 



edk) = I 1 

—ivkR 



Oik^), 



where ly is given in Eq. (|12p and vl = 1.99 x 10"' m/s. To 
lowest order in eL{k) points along the CNT axis Cy, 
as expected for a longitudinal mode. Finally, the radial 
breathing mode corresponds to 



UJB{k) 



iK + p 




O(fc'), 



(27) 



esik) = ( -ivkR | +0{k^). 

This mode has an energy gap, Hlub — 14: meV for R = 
1 nm, scaling a,s lub ex. R~^. The quoted results for the 
velocities vt,l and the frequency ujb, first obtained in 
Rcf . m, follow from Eq. pT|) , and are in accordance with 
ab-initio calculations.— 

For angular momentum £ = ±1, we recover the correct 
dispersion relation of the important flexural (J — F) 
modes.— They arc degenerate and correspond to 



ujp{k) 
eF,t=±{k) 



1 

71 



1 

TkR{l 
T^ (l- 



V2vlR' 

4 \ 



(9+Qiy)k^Ft- 
4 



(28) 



0{k^). 



J 



Note that for |fc|i? < \/2vt/vl, and thus for all wave- 
lengths of interest here, the flexural phonons arc the 
lowest-lying modes available. 

Next we observe that for £ 7^ 0, longitudinal modes 
acquire a gap, a'L(fc ~Q,£) — vt\£\/R, and "breathing" 
modes have an even larger gap than Eq. (P7)) . cjb(0, £) = 
vT+7^a;B(0, 0). Since we focus on low-energy acoustic 
modes, these gapped modes are irrelevant and will not 
be studied further. Moreover, the diagonalization of the 
elastic matrix (jl9p shows that for any £ flexural modes 
remain gaplcss. However, for \£\ > 1, curvature effects 
(see App. [A| will open gaps for these modes as well.— 
For i? ^ 1 nm, such gaps are comparable in magnitude 
(or slightly smaller than) the frequency of the breath- 
ing mode ([27]) .— Since ph-ph interaction effects become 
more and more pronounced with decreasing radius R (see 
below), the most interesting application range of our the- 
ory is 4 A< i? < 1 nm, where most phonon modes have 
rather large gaps but a continuum elasticity approach 
is still reliable. Ignoring gapped modes is then a good 
approximation over a wide temperature regime and in 
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our low-energy approach we need to retain only gapless 
modes, i.e. the T mode (P5|) . the L mode (pS)) . and the 
two degenerate flexural F modes (f28|) . The gaplcssncss 
of the £ = ±1 flexural modes is robust against curva- 
ture effects and protected by rotational symmetry. Note 
that the resulting theory is only valid on energy scales 
below those gaps; for R w 0.5 nm, this is justified up to 
temperatures of order of 50 K. 

From now on, the sums over {J,£) will then only run 
over (T, 0),(L,0), and (F, ±). It is remarkable that for 
all these phonon modes, elastic continuum theory is able 
to yield accurate dispersion relations which are in good 
agreement with elaborate force-constant^^i^i^i^^ and ab- 
initio calculations.— Since the breathing mode ([27|) may 
be of interest for future thermal expansion calculations, 
we specify the corresponding three-phonon matrix ele- 
ments in App. [B1 but for the main part of the paper we 
will neglect this mode. 



III. PHONON-PHONON INTERACTION 
PROCESSES 

In this section, we evaluate the three- and four-phonon 
scattering amplitudes following from Eqs. ^5)) and ([T7]), 
respectively. They are obtained by inserting the normal- 
mode expansion for the displacement field into the 
definition of the strain tensor, see Eqs. (O and We 
will always keep the lowest nontrivial order in <C 1, 
but also specify the next order when cancellation effects 
are anticipated for the leading order. It is then straight- 
forward to obtain the full normal-mode representation of 
the nonlinear strain tensor. The result can be found in 
explicit form in Appendix [B] 



dcr in |fc^|i? < 1, 



A 



LLL 



= --{l-v)kik2kz[2K 

+ (1 - + (1 + ^.2)^^] 



(30) 



Altt = -2([2M + e2(l-!^)]fcifc2fc3 (31) 
- v^ikfj, 

-iSe^-esfJ-i^ + '^)kik2k3, (32) 

-^-^^Si2^^i3kik2k3{k2 ~ k3)R, (33) 



AL,Fe2,Fe3 

At.F12,F13 



with V in Eq. (fT^ . The matrix elements related to 
the breathing mode can be found in App. [BJ Sym- 
metry under phonon exchange is taken into account in 
the expressions ((5n)) - (P5|) . and the {Ji,J2,Ji) summa- 
tion in Eq. ^ runs only over (LLL), (LTT), (LFF) 
and (TFF) , while all other matrix elements vanish iden- 
tically. In particular, there is no amplitude for LLT 
processes^ nor for the scattering of three twist modes, 
TTT . Moreover, all amplitudes involving an odd number 
of flexural phonons vanish by angular momentum conser- 
vation. We also observe that the anharmonic third-order 
couplings ^1 and ^2 do not introduce new physics, but 
only rcnormalize parameter values of coupling terms gen- 
erated already by geometric nonlincaritics. In fact, the 
leading contributions to phonon decay rates turn out to 
be completely independent of such anharmonic couplings, 
as we will show in Sec. |lVl see Eqs. (IMl), (ES) and (p7|) 
below. 



B. Four-phonon processes and flexural phonon 
interaction 



A. Three-phonon processes 

The normal- mode representation of the strain tensor in 
App. IB] allows us to write iJ*-^-* from Eq. pB]) in the form 
of a standard three-phonon interaction Hamiltonia n^^i^^ 
(note again that Y.k = IZo If)' 



1 



J2 Y.'^JM(.ki,k2,k3) (29) 



J1J2J3 kik2 

X Uj^{ki)uj^{k2)uj3{k3), 



Next we turn to four-phonon interactions. A similar 
result as for three-phonon interactions, see Eq. can 
be derived using the strain tensor given in App.lBl Since 
FFF matrix elements vanish, quartic terms are crucial 
in the case of flexural modes, and we shall only discuss 
these four-phonon matrix elements in what follows. It 
is nevertheless straightforward (if tedious) to study also 
other four-phonon matrix elements based on the expres- 
sions given in App. [B] 

Since the flexural mode is the lowest-lying phonon 
branch, FFFF processes provide the only possibility for 
its decay at T = 0. It turns out that the relevant coupling 
strength for such processes is parametrized by 



where the £ summation is implicit when J ^ F, i.e. J 
stands for both the phonon mode index and the angu- 
lar momentum £. Due to momentum conservation k^ = 
— fci — k2, and uj{k,£) has been defined in Eq. (P^ . Af- 
ter some algebra, we obtain the following non-vanishing 
three-phonon amplitudes Aj-^j^j^ (fci, ^2,^:3) to leading or- 



K 



2ttR 



m \ ' 
PmJ 



(34) 



where K = fi + \ and m is given in Eq. ([28|) . Note that 
g (X 1/R^, and thus flexural phonon interactions become 
stronger for thinner CNTs. 
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After some algebra we find 



AFei,Fe2,Fe3,Fe4.{ki, k2, k^, k^) 



2ttR ^ 

felfe2fe3 {^} 

X UFiki,£i)uF{k2,i2)uF{k3,e3)uFik4:, 4), 

where ^4 = — (/ci + fc2 + fca), angular momentum conser- 
vation implies the condition ii + £2 + + £4 = 0, and 

fcifc2fc3fc4 



Fli,Fl2,Fl3,Fti 



K 



Now for all {£i\ combinations with X]i=i = and = 
±1, one finds X]i<j ^j^j = ~2 and J|,- = i. This allows 
us to carry out the £ summation, and gives 

H'^'^^ = g{PM/mf kik2k3ki (35) 

X UF{ki,+)uFik2,+)uFik3,-)uFik4, -). 

Note that Eq. ([55]) is determined by geometric nonlinear- 
ities alone, i.e. by the contribution of the nonlinear part 
of the strain tensor in U2 [u] . The anharmonic third- and 
fourth-order couplings (^1,2 and ki,2,3, respectively) also 
give rise to Affff contributions, which however con- 
tain higher powers in \ki\R <^ 1. Such anharmonic four- 
phonon processes arc therefore parametrically smaller 
than the geometric nonlinearity (|35p . and can be ne- 
glected in a low-energy approach. 

In coordinate space, Eq. (j35p corresponds to a local 
four-phonon interaction. To see this, we represent the 
momentum conservation constraint in Eq. psp as 

S{ki +k2 + k3 ■ ^ '^^ 

and then arrive at 



-i(ki+k2+k3+ki)y 



2ti 



H = dy 



P 



Qy2 Qy2 



+ giPAi/m) 



Pm 4to2 

2 

2 I uu_uU 

dy dy 



(36) 



where pM = '^t^Rpm is the effective linear mass den- 
sity, and the (non-Hermitian) coordinate-space flexural 
displacement operator is defined as 

u{y) = -±=Y.e^''yuF{k,+) (37) 



2pMUJF{k) 



Jky 



(aj,(fc,+) +a],(-fc,-)) , 



hpMUJFik) 



with the canonically conjugate momentum field operator 
p{y) = \/2^^e*'=^pj.(fc,-) = 

k k 

e^'^y (^aj,{k,-)-aU-k,+)). 



X e 



Since g > 0, the interaction among flexural phonons is 
repulsive. Therefore phonon localization and two-phonon 
bound states^ are not expected to occur. 

Remarkably, as we show in detail below, it turns out 
that the finite-temperature decay rate TF^k) for a fiexu- 
ral phonon diverges when the interaction (j35p is treated 
perturbatively. A related breakdown of perturbation the- 
ory for phonon decay rates has also been reported by 
PerrinS in a study of optical phonons in molecular crys- 
tals. In that case, the singularity could be traced to the 
flatness of the dispersion relation. A similar situation 
occurs for the magnon decay problem in ID spin chains, 
where the analogous perturbation theory also predicts 
a finite and momentum-independent T = decay rate 
above a certain threshold, while the correct (nonpertur- 
bativc) result vanishes at the thresholdsi^ In our case, 
the divergence arises due to the conspiracy of the almost 
flat dispersion relation, uJF(k) = hk'^ /2m, with the low 
dimensionality (ID). This implies a macroscopic phonon 
generation in the noninteracting case for finite T . For 
a system of length £, the total number of £ = ± ficxu- 
ral phonons follows with the Bose-Einstein distribution 
function (/? = l/fcsT), 



n{uj) = 



1 



(38) 



as = 2£ ri(a'i?(fc)). As a result, the ID phonon 
density p ~ N/C ~ 2mkBTC/TTh^ diverges in the ther- 
modynamic limit C —I- 00 any finite temperature T. 
This situation therefore calls from the outset for a non- 
perturbative treatment of the interaction (j35p . In view 
of the divergent noninteracting phonon density p, we ex- 
pect that mean-field theory is able to properly handle 
the regularizing effect of the interaction despite the low 
dimensionality, at least in a semi-quantitative fashion. 
We thus employ mean-field theory to compute the ID 
flexural phonon density p{T), and then use this result in 
Sec. llVl for the decay rate calculations. 

Taking fike = {a}p{k,£)ap{k,£)) as the only non- 
vanishing mean- field parameters in Eq. (j35p . the mean- 
field Hamiltonian for the flexural modes is given (up to 
irrelevant constants) by 



kl=± 



[nwF{k)+^gp]o^p{k,£)aF{k,£). (39) 



where p = Ylki^kt- In order to derive Eq. ([39]) . 
we disregard all terms involving an unequal number of 
creation (a) and annihilation (a^) operators, see also 
Ref. [5^. Moreover, for all nonvanishing contributions 
to the mean-field approximation of Eq. (j35p one finds 
sgn(fcifc2fc3fc4) = 1. The resulting self-consistency equa- 
tion is then p = 2 n{LOF{k) -f 4g/9). The momentum 
integral can be carried out and yields 



P 



ksT* 
43 



{T/T*)Y{T/T* 



(40) 



The temperature dependence of p shows universal scaling 
with X = T/T*, where we introduce the temperature 
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scale 



T* 



(41) 



The dimensionless scaling function Y{x) is determined 
by the self-consistency condition 



Y = Lii/2 (e-^) , 



(42) 



with the polylogarithmiai Lis(z) = J2i=ij Equa 



tion p2|) can be analytically solved in the limits x <^ 1 
and x 3> 1, and allows for numerical evaluation in be- 
tween those limits. In particular, we find Y(x ^ 1) ~ 
— ^ln(7rx) and Y{x ^ 1) ~ where we exploit the 

relation^ limy_»o Lii/2(e^^) = ^J-k/Y . In accordance 
with our above discussion, we therefore find that in the 
noninteracting (T* = 0) case, p diverges for any finite 
T . However, once interactions are present, a finite flexu- 
ral phonon density p emerges, which for T ^ T* can be 
written as 



p{T > T*) = MI(T/r*)2/3. 

4g 



(43) 



Using the parameters in Eq. pT|) . the scale T* in Eq. (gT 
is estimated as 



3.7 X 10-9 K 
(i?[nm])7 ■ 



(44) 



Even for the thinnest possible CNTs (where R w 0.3 nm), 
this puts T* deep into the sub-milli-Kelvin regime. As- 
suming T ^ T* from now on, we take p as given in 
Eq. Within mean- field theory, see Eq. ([35]), non- 

perturbative effects of the interaction ([35]) thus lead to 
the appearance of a dynamical gap tUp = 4gp for flexu- 
ral phonons. We effectively arrive at a modified flexural 
dispersion relation^ 



a;^(fc)=a;p + — , 



(45) 



characterized by the temperature-dependent gap 

fy^p{T) = {T/T*f'^kBT*. (46) 

From now on, we take Eq. (|45p for the dispersion rela- 
tion of flexural phonons. Since T ^ T* , we also observe 
that the gap is always thermally smeared, kgT ^ hiOp. 
Nevertheless, it is crucial when discussing the decay rate 
for a flexural phonon. 



IV. DECAY RATE AND QUALITY FACTOR OF 
ACOUSTIC PHONON MODES 



In this section, we study the decay rate of a phonon 
excitation with longitudinal momentum p = hk > and 



mode index J = L,T or F. We compute the finite- 
temperature decay rate Tj{k,£) from lowest-order per- 
turbation theory in the relevant nonlinearity. At T = 0, 
this corresponds to the standard Fermi's golden rule re- 
sult. 

A. Self-energy calculation 

To access the finite-T case, we will first write down the 
respective imaginary-time self-energy T.j(T,k,£), where 
< T < hf3 denotes imaginary time. The Matsubara 
Green's function is defined via 

(uj{nn,k,£)uji-n'„,~k\~e)) = 

-27rhpSa„,a'Jik - fc')Gj(f^„, fc, ^), 

where the r2„ = 2-Kn/h[3 (integer n) are bosonic Mat- 
subara frequencies and uj{t, k,£) is defined in Eq. (j24[) . 
Employing Eqs. ((22)) and (|24l) . the noninteracting Green's 
function is 



G^°\nn,k,£) 



'fl/PM 



nl+u;jik,£) pM 



g{n^,uj{k,£)). 



The function ^y(f2„,u;i) has the time-representation 



(47) 



(48) 

with the Bose function ([55)) . The full retarded Green's 
function Gj(uj,k,£) follows after analytic continuation, 
ifln —^L^J + zO^, with the Dyson equation, 

Gj\co,k,£) = {G^°^)j\u;,k,£)~^j{Lu,k,£), (49) 

leading to the on-shell [uj = Luj{k,£)] rate 

h 



Tjik,£) 



Pm^ 



-ImSjluJ, k, 



(50) 



The relevant quantity needed to estimate the decay rate 
is therefore the self-energy T,j{uj,k,£), whose imaginary- 
time version is Sj(r2„, fc, £). 

In the self-energy calculation for the various modes 
shown below, we will encounter integrals of the type (in- 
teger r > 1) 

/,(r!„;^i,...,a;.) = / dTe'''-^Y[g{T,uj,). (51) 
Jo 

Empl oying Eqs. ([i5|) and ([55]) . we find for r = 2 (see also 
Refs. \mMi 



/2(rJ,i, wi,a;2) = E 



l-f n(^iwi) -I- ^(^2^2) 



Similarly, for r = 3, we obtain^ 



iuJiLJ2 i^n + ^It^l + ^2^2 



(52) 
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E 



n(^iWi)n(^2W2)n-(6w3) 



r 



(53) 



Let us then proceed with the discussion of the different 
phonon modes, starting with J = L. 



where I2 is given in Eq. ([5^ . For kR <C 1, after analytic 
continuation, we obtain the rate from Eq. (|50p . Identify- 
ing u! = VLk yields with Eq. ([55]) the result 



B. Longitudinal stretch mode 



The dominant contributions to the decay rate for a lon- 
gitudinal phonon come from the relevant non-vanishing 
three-phonon matrix elements, namely L L + L in 
Eq. dSOl), L ^ T + T in Eq. and L F + F in 

Eq. ([32]). The amplitude for the i ^ i + T process 
vanishes, and such decay channel could only be possible 
via higher-order processes involving the virtual excita- 
tion of flexural modes. While one could compute the 
corresponding contribution, we expect that it is negligi- 
ble against the rate found below. We also anticipate that 
the contribution of the process L ^ T + T is sublead- 
ing with respect to that L ^ F + F, as dimensional 
arguments at T = suggest and the explicit finite-T 
calculation shows. Therefore, we also neglect this decay 
channel. Finally, while the process L ^ L -I- L is in prin- 
ciple kincmatically allowed for a strictly linear dispersion 
relation, energy conservation cannot be satisfied as soon 
as one takes into account the 0{k'^) corrections to ujL{k). 
Thus, this decay channel can also be safely omitted. 

The only remaining possibility is then the process 
L F + F, where the two flexural phonons carry oppo- 
site angular momentum. Energy conservation then poses 
no problem as long as Wp ^ wl/c, see Eq. (^5)) . For clarity, 
we now focus on this case, where the channel L ^ F + F 
provides the dominant decay mechanism for a L phonon. 
The lowest order in perturbation theory generating a fi- 
nite decay rate comes from the "bubble" diagram (i.e. the 
second order). 



Sl(t, k) 



^\ALFF{k,qi,q2)\'^giT,uJi)g{T, UJ2): 



where qi^2 = T9 + k/2 and u;i_2 — ^'i'i2/2"T'i and 
the amphtude Alff in Eq- ([32]) is evaluated, say, for 
I2 = —£3 = 1. The two Green's hmctions correspond to 
flexural phonons. Using ii{\ + v) — pMv\l2, we then find 



EL(ri„, k) 



{hvLkf 



dq UJiUJ2l2{^n]^l-,^2)-, 



rL(fc) 



dq 



{[n{uji) + n{uj2) + l] 



uji — UJ2) + 2[n{uj2) — n{uji)]5{uj — ui + CJ2)}- 
We now need to resolve the (S-functions represent- 
ing energy conservation. The first term yields q = 

1 



± 



— while the second leads to q ^ 

We then collect terms, keeping only the leading order 
in kR <Si 1 and ksT ^ hv^/R - otherwise the thermal 
scale ksT would exceed the smallest gap of the discarded 
phonon modes, and we would need to account for the ef- 
fects of the mean-field gap tUp. We find 



rL(fc) 




coth(fi/3wifc/4) (54) 



smh{hPvLk/2) 



Since this rate does not depend on the anharmonic cou- 
plings (^1,2 and ^1^2,3), we find the remarkable result that 
the dominant decay rate for the longitudinal phonon is 
solely determined by the Lame coefficients (or, cquiva- 
lently, by the sound velocities). For T = and k ^ 0, 
the rate becomes universal, i.e. completely independent 
of material parameters. This is due to the fact that in 
our elastic model the curvature of the flexural branch 
~ l/m is proportional to the longitudinal sound velocity 
vl- Moreover, the T = rate is cx ^/k, i.e. the nonlinear 
damping effects become important at long wavelengths. 

For an estimate of the quality factor for longitudinal 
modes, we now put k = tt/C in Ql = uj ^{k) /T ^{k) , with 
CNT length C. This yields 



Ql{T,C) 



coth 



\ACkBT 



\2CkBT ) 



(55) 
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with the T = result 



Ql(0,/:) = (3.80 X lO*') X {R[nm\f y^R/C, (56) 



([TT]) . For typical parameters, say, 
= 500 nm, this gives Ql{T ^ Q) ^ 



where we used Eq. 
R ~ 0.5 nm and C 
1.5 X 10''. 

This number for the quality factor implies a sur- 
prisingly strong damping effect due to phonon-phonon 
interactions, and is similar to what is observed 
experimentally^^ Note that this value yet has to be un- 
derstood as upper bound since there might be other de- 
cay mechanisms. Since LOp{T = 0) = 0, inclusion of the 
mean-field gap Wp in the above derivation does not affect 
the estimate ([5^ . Furthermore, with Ql{Q,C) oc 
we observe that for sufficiently long CNTs and low T, the 
anharmonic decay is always important, in agreement with 
the conclusions of Ref. [11, see our discussion in Sec. H] 
The temperature dependence of Ql is controlled by the 



J 



ratio of the longitudinal confinement energy scale, Hvl/ C, 
and the thermal energy, fc^T. For ksT ^ Kvl/ C, we find 



Ql{T,C) 



^flVL ^ R"^ 



CZ/2T' 



(57) 



We observe that the damping effects get stronger for thin- 
ner CNTs. 



C. Twist mode 

Next we turn to the twist (J = T) mode, where the 
vanishing of the T ^ T + T amplitude and the absence of 
the (kinematically forbidden) T T + L channel imply 
that the T F+F decay will dominate. The calculation 
proceeds in the same way as for the L phonon, and in the 
low-energy limit, ksT ^ Hvt/R and ujp <C vrk, we find 



Trik) 



h VT ,/^ 2^/\kRfr 



(58) 



+2^/'^[{vL/vT)kR]^'^/^e ■^v^.r^Rk^T gj^h 



(hvTk ^ 
2k^^ 



We note that the T = rate is cx A:^/^. 



r 



Putting k = n/C, the quality factor is then 

Qt(T,/:) = -^. (59) 

This gives for T = the estimate 

Qt(T = 0,/:) ~ (5.74 X 10^) X {R[iiia]f y/C/R, (60) 

showing that damping of the twist mode, which is en- 
ergetically below the L mode, vl/vt — 0.62, is much 
weaker, in accordance with Ref. [4^ In the high- 
temperature hmit, we find Qt oc i?^/^/(r£^/^). 



in Figure [TJ Note that we have already taken into ac- 
count the self-consistent tadpole (first-order) diagram by 
using mean field theory, leading to the dispersion rela- 
tion (|45p . The imaginary-time self-energy corresponding 




FIG. 1: Self-energy diagram contributing to the decay rate 
Tpik) in leading order ("fishbone diagram"). Solid lines de- 
note a flexural phonon propagator, taken from mean-field the- 
ory. 



D. Flexural mode 



Next we discuss the decay rate Fi?(fc) for the flexural 
mode. We put £ ~ +1 (the rate for ^ = — 1 is identical), 
and describe the perturbative calculation of F^^ (k) . The 
leading term involves the decay of the phonon into three 
flexural phonons, F ^ F + F + F, see Eq. (|35|) . corre- 
sponding to the "fishbone" diagram for the self-energy 



to the fishbone diagram involves three Green's functions, 
with momenta qi and frequencies Ui = ujp + hqf/{2m), 
where z = 1, 2, 3. Employing Eq. (|53p and the interaction 
strength g in Eq. p4p . we then obtain 
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S^(17„, k) = [ dq,dq,dq,^MsLs(q, + q, + q, - k) (61) 

47r^m* J UJ1UJ2UJ3 



, , , . i^n + Cl^l + £.2^2 + 6^3 + 6t^2 + ^3^3) ' 



E 



The (5-function represents momentum conservation. After analytic continuation, using the relation 2n{±uj) = 
± coth(fi,/3w/2) — 1, this yields the rate 

rF(fc) = jT— — ^ y_] / dqidq2dq3d{qi+q2+q3~k) (62) 



<5(a. + eicci + 6^2 + 6^3) n ( -[coth(;i/3c^,/2) - 



i=l 
I 



We now express the (5-functions as 



27r 

dt 
2^' 



^+Ei?.'^i)* 



which decouples the qi integrals and allows to perform 
the summations. Setting uj = ^^(fc), the on-shell rate 
reads 

r-(^) - 77^.1^ / dtdye-^^'^y'-^^G^it,y), (63) 

^^[1 + 71(0")] J 

with the correlation function 

dq hq^ cos{qy) 



G{t,y) 



(64) 



m 27r 2pMiOF{q) 

X (cos[a;F(9)^] coth[?i/3a;j?(q)/2] — i sin[ti>i?((7)t]) 

m 



= ^([a,^.t](i,y)[a,^,](o,o))o. 



Here u{t, y) is the Heisenberg representation of the flex- 
ural displacement operator u{y), see Eq. p7p . and the 
noninteracting average (■ • • )o is taken with respect to 



fl"MF, see Eq. The product of three Green's func- 

tions in Eq. ([55]) reflects the structure of the fishbone 
diagram in Fig.[TJ Using the dispersion relation (1451) and 



in Eq. (|46)) . we observe that the g-integral for G{t,y) 
is regular. 

The rate T p for the decay of the mode with wavevector 
k = tt/C then depends only on the two dimensionless 
quantities 



C 



T 



where we define the lengthscale 



C* 



(65) 



(66) 



\/2mfc^T* 8mg 
Using Eq. (jH]), we obtain the estimate 

C*[^lTn] ~ 534 (i^[nm])^ 

which gives C* — 33.4 /im for R = 0.5 nm. By rescaling 
all lengths in units of £* and all frequencies (or inverse 
times) in units of ksT* /h, Eqs. and ([M]) imply a 
universal result, where the dependence on material pa- 
rameters only enters via T* and C* . Using Eqs. P5|) and 
after some algebra, we obtain 



kuT* 



1^2/3 y2 



dydt exp 



r -y2/3 



(67) 




2/3 27rti; 



Aj,' cos 



v2/3 

XrJ y 



[cos(u'<) coth(u)/2XT) — isin(zi;t)] 



For T = 0, Eq. (p7|) can be solved in closed form and (In fact, this result easily follows also from Eq. (p^ .) Re- 
gives markably, this rate does not depend on the length {Xc) of 

rHx,.x, = o, = ^£i:. (68, 



12 




R[nm] 

FIG. 2: Quality factor for the three low-energy modes at T = 
and £ — 500 nm as a function of tube radius R. The 
solid line depicts Qf, the dashed line Ql and the dotted line 
10-^Qt. 



the CNT, i.e., it is independent of phonon momentum k. 
Despite the smallness of T* , see the estimate in Eq. (|44)) , 
this predicts a surprisingly large damping effect due to 
phonon interactions. Estimating the zero-temperature 
quahty factor as above, we find 

QpiT = 0,/:) ~ (1.98 X 10^^) X {R [iiiii]f{R/Cf. (69) 



Taking C = 500 nm and R = 0.5 nm, this gives Qf ~ 
3 X 10"^. 

The i?-dependence of the zero-temperature quality fac- 
tors for the various modes is summarized in Fig. [2l In 
the range of radii considered, it is much stronger for the 
flexural mode than for the longitudinal and twist modes, 
which have an approximately similar dependence. Note 
however that Qt is in fact three orders of magnitude 
larger than Ql- 

However, as discussed in Sec. IIII Bl due to the small- 
ness of T* , the zero-temperature limit is actually not ac- 
cessible experimentally, and one is always in the regime 
T ^ T* . It is then interesting to evaluate the rate in 

1 /3 

the limit Xt ^ 1 but with Xj! Xc ^ 1, which corre- 
sponds to a regime in which the kinetic energy is much 
larger than the gap. For sufficiently short CNTs (or if one 
evaluates the rate at a larger wavevector than 7r/£) the 
two inequalities can be satisfied simultaneously. In that 
case, we find from Eq. ^ that hVp/kBT* oc X^X^. 
The resulting power-law behavior for the tempera- 
ture dependence of the flexural phonon decay rate is a 
prediction that should be observable with state-of-the art 
experiments. 



V. ELECTRON-PHONON COUPLING 

In this section, we consider the effects of e-ph couplings 
in metallic single-wall CNTs. In contrast to the semi- 
conducting case studied before, the coupling of phonons 
to electrons may contribute another decay channel be- 
yond ph-ph interactions. Within lowest order perturba- 
tion theory, the two mechanisms are additive for the de- 
cay rates, and hence also for the inverse quality factors. 
Here, we only take into account the deformation poten- 
tial as a source for e-ph coupling, as this coupling has a 
very significant strengthi^i^i The phonons, described by 
the strain tensor Uij^ create the deformation potential^^ 

Ve-M = ^oTr [^""(r) + ^"''"(r)] , (70) 

where Vb ~ 20 eV (see Ref. |33| ) . We will consider spinless 
electrons at a single Fermi point and eventually multiply 
by a factor 4 the final result for the rate, to take into 
account the electron spin and KK' degeneracies. More- 
over, we only take into acccount electrons in the lowest 
transverse subband, with zero angular momentum. This 
is justified since the higher angular momentum states are 
separated by a large energy gap. Then, from the normal- 
mode representation in Appendix[Bl we observe that only 
the L mode couples via the first term (u''") in Eq. ([70]) 
(the coupling to the F mode requires at least one electron 
in a higher angular momentum state) , while the F and T 
phonons couple only via the second (u'^''") term, where 
two phonons are involved. The electronic contribution 
to the decay rate of the T phonon is then expected to 
be weak, and we will focus on the decay of the L and F 
modes. 

Let us start with the L phonon with momentum hk, 
where the relevant lowest-order self-energy diagram due 
to the first term in Ve^ph corresponds to the electron-hole 
bubble. Its imaginary part, responsible for the phonon 
decay rate, is (x [5{uj~vpk) — 5{uj-\-VFk)], where to = VLk. 
(Note that we use a linearized dispersion for electrons.) 
Since the Fermi velocity is vp ~ 10^ m/sec ^ vl, this 
condition can never be met, and only higher-order con- 
tributions can possibly lead to a phonon decay. This sug- 
gests that the decay rate of the L phonon due to Eq. ([7D1) 
is very small, and probably negligible against the ph-ph 
mechanism. 

For the case of a F phonon, the lowest-order contri- 
bution comes from the second term in Eq. ([7CT)) . lead- 
ing to a "fishbone" diagram with two electron lines and 
a phonon line. The corresponding imaginary-time self- 
energy is given by 

(71) 

with loq = ujp + h{k — qY /2'm and the ID electron polar- 
ization functions^ 

ne(r,9) = ^5^n(Ca.e)e«"=^ 
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where We = i'fI'?!- Following the same steps as in Sec. lIVl 
we then find for the T = decay rate 

1 

which yields the energy conservation condition 

vlR 

For q 0, this condition can only be met if kR > 
vf/{V^ vl), i.e. only for short-wavelength phonons. 
Thus, for the long-wavelength phonons of interest here, 
the energy mismatch between electron-hole pair excita- 
tions and the phonon modes implies that again only 
higher-order terms can possibly generate a finite decay 
rate. 

The above discussion therefore suggests that e-ph cou- 
plings via the deformation potential do not lead to sig- 
nificant decay rates of the gapless L and F modes. Their 
decay should then be dominated by the ph-ph interaction 
processes as described in Sec. IIVI 



to important consequences and allows to compute a well- 
defined decay rate for flexural phonons in low-order per- 
turbation theory. 

Using this approach, we have determined the decay 
rate and the quality factor for all long-wavelength gap- 
less phonons in carbon nanotubcs. We have also shown 
that elcctron-phonon interactions arc ineffective in re- 
laxing those modes due to a mismatch in energy scales. 
The reported quality factors (Q) are remarkably small, 
especially for thin CNTs, pointing to important phonon 
damping effects. Phonon-phonon interactions in CNTs 
are therefore significant and quite strong. Moreover, the 
values we have found are rather close to the typical Q re- 
ported in recent experimentsii^ Our predictions represent 
intrinsic upper bounds for Q. Such upper bounds can be 
valuable in assessing the predictions of approximate the- 
ories, or when interpreting experimental data in terms of 
phonon damping. We hope that our work will motivate 
further experimental and theoretical studies along this 
line. 
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VI. CONCLUSIONS 

In this paper, we have formulated and studied a gen- 
eral analytical theory of phonon-phonon interactions for 
low-energy long-wavelength acoustic phonons in carbon 
nanotubcs. The continuum elasticity approach employed 
here reproduces the known dispersion relations of all gap- 
less modes, including the flexural mode, ujp{k) = hk^/2m 
with "effective mass" m. We have then included the most 
general cubic and quartic elastic nonlincarities allowed 
by symmetry. Remarkably, the relevant phonon-phonon 
scattering processes giving the dominant contributions 
to the decay rates are already found from the geometric 
nonlincarities (i.e. using the nonlinear strain tensor in 
a lowest-order expansion of the clastic energy density), 
and for a quantitative discussion of the decay rates, only 
the knowledge of the well-known Lame coefficients (or 
equivalently of the sound velocities) is necessary. 

We have provided a complete classification of all pos- 
sible three-phonon processes involving gapless modes, 
along with the respective coupling constants. At the level 
of four-phonon processes we have focussed on the flexu- 
ral modes, where a peculiarity is encountered, since the 
four-phonon processes lead to a singular behavior of the 
finite-temperature decay rate. The physical reason for 
this singularity is the proliferation of phonons at finite 
temperature, implying a divergent ID density of flexural 
modes. Interactions effectively regularize this divergence 
and lead to a finite density. We have employed mean- 
field theory to quantitatively describe this mechanism, 
and found a dynamical temperature-dependent gap fkUp 
for fiexural phonons. While this gap is in practice al- 
ways below the thermal scale ksT, it nevertheless leads 
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APPENDIX A: CURVATURE EFFECTS 

In this appendix, we briefly discuss how to include cur- 
vature effects in the nonlinear elastic energy. To that end, 
we consider the metric tensor for the cylindrical surface, 
whose components j ~ x, y — 1, 2), 

n—x^y,z 

are expressed in terms of the ordinary scalar product of 
tangent vectors 

E9x(r) 
^n,i (^^)*^n (^) — 7^ : 
OXi 

n 

with X given in Eq. The nonlinear strain tensor 

then follows equivalently from 2uy = gij — g^^^ , with the 

metric tensor g^^^ = 6ij of the undeformed cylinder. The 
local unit normal vector is N = (tj. x tj,)/|t3; x tj,|, and 
the mean local curvature Q of the cylinder is defined as 

where g^^ is the inverse of gij , and the second fundamen- 
tal form of the surface is 

bij = ^ NndjtnJ. 

n—x,y,z 
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To lowest order in the displacement fields u„, the mean 
curvature is then given bj*^ 

where JIq = —1/2R is the curvature of the undeformed 
cylinder. 

Curvature leads to an additional energy cost due to hy- 
bridization effects. One can model this in a phenomcno- 
logical way by adding an additional term 

Z^curvN = K-i^ - ^of 

to the elastic energy density, where k is a proportion- 
ality constant. Such effects turn out to be small un- 



less one deals with ultrathin CNTs, but they provide 
gaps to flexural modes with angular momentum \£\ > 1. 
For i? < 1 nm. these gaps are typically comparable in 
magnitude^ to the breathing mode energy in Eq. ^7} . 



APPENDIX B: NORMAL MODE 
REPRESENTATION OF THE STRAIN TENSOR 

In this appendix, we provide the explicit form of the 
strain tensor expressed in terms of the normal mode dis- 
placement field operators uj{k,£), see Eq. ([24|) . We keep 
all ^ = modes {L,T,B), and the gapless flexural (F) 
modes with £ = ±1. For the linear part of the strain 
tensor, see Eq. (O, we obtain from Eq. ([23]) the result 



it""(r) 



V2^ 
ik^R 



,iky 



-iuk , ,,s 
^k ) 



ik/2\ 
zfc/2 )"^(^') 



1/i? 
vPR 



UB{k) 



V2 



2^^ I {l + v)kR ~i 1''^^^^ 



while the nonlinear part ^ reads 
u"""(r) 



2ttR f-^ 1 V -kk^U 



ui^(ki)ui^(ki) 



^fci(fc2 - kx) \ /, \ n \ 

Pi(fc2-fci) ' )^L{h)ur(k,) 

+ I ^ [uT{kl)uT{k2) + UB{ki)uBik2)] 



;(fcl - k2) 



-ivki/R 



+ 



V2 



E 



1/2^2^1 (l+Z^)fcl/c2/2 \ . „ . 

(l + .)fclfc2/2 £2(fc2+-fcl)fclfc2i? j"^^'^^)"^^^^'^^) 



vkl -l2k2/2R 

i2k2/2R -fclfc2 



UT{ki)uF{k2,£2) + 



iv£2kl ik2/2R \ 
ik2/2R ii2kik2 ) 



UB{kl)uF{k2,£2) 



kik 



\^ „^x{i,+i2)/R 1 2 

A I ieiki+e2k2+H(i+i2){ki+k2)]R (^ - p p \ 

£ie2 \ 2 <^i<^2; 



y.uF{ki,li)uF{k2,t2) >■ 



For completeness, we also list the long-wavelength form mode: 
of the three-phonon amplitudes involving the breathing 

Abbb 
Abtt 



2i?3 

2i?3 ' 

i?2 



^LBB = -i[-/C + M + 3(l-!^Ki + (l-3i.)6] 



Alls = -'^^K + 3{1 ~ ,.')(^ + {2,, - 1 - 3i^')U 

jSK - fi)k2k3 

Arff = —Olo.-t, tt: • 
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